clear all
clc
load closure.mat
% 设多元线性回归方程 y=a1*x1+a2*x*+a3*x3+a4*x4+a5*x5+c
x1 = data(:,1)
x2 = data(:,2)
x3 = data(:,3)
x4 = data(:,4)
x5 = data(:,5)
y = data(:,6)

x = [ones(size(x1)),x1,x2,x3,x4,x5]
[b,bint,r,rint,stats] = regress(y,x)
disp('多元回归方程常数项')
disp(b(1,1))
disp('多元回归方程系数项')
disp(b([2:6],1))
disp('***********************************************')
disp('是否绘制残差图像？输入0否，1是： ')
temp = input('')
if temp == 1
    rcoplot(r,rint)
else
end
